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The single channel transport properties of a quantum wire network in an ambient electromagnetic 
field are investigated. The network is made of three semi-infinite external leads connected to a ring 
which is crossed by a magnetic flux. Two different kinds of bulk dynamics, associated with the 
Schrodinger and Dirac equations, both minimally coupled to the external electromagnetic field, 
are analyzed. Connecting the external leads to heat reservoirs with different temperatures and/or 
chemical potentials, the system is driven away from equilibrium. We derive the conductance and 
the noise of this configuration, describing in details the impact of the magnetic flux. In the case of 
weak coupling between the ring and the reservoirs, a resonant tunneling effect is observed. We also 
discover a new power law behaviour of the pure thermal noise which differs from the usual linear 
Johnson-Nyquist law when the magnetic flux is nonzero. 



I. INTRODUCTION 

Quantum wire networks are expected to play a funda- 
mental role in future nano-electronics. For this reason 
the transport properties of such networks attract much 
attention P-[35|. Since a quantum wire network is com- 
posed from junctions and extremely thin (of the order 
of few nanometers) leads, it is natural to represent it 
by a graph, as shown in Fig. [ij where the vertices Vj 
and edges Ei model the junctions and the wires respec- 
tively. In this way one is led to the study of electron 
transport along graphs. This problem has been exten- 
sively studied in the last three decades (see for instance 
|36)37j and references therein) with a recently growing 
interest in the universal features of the above systems. 
Despite great progress in this direction, the role and the 
physical implications of the different possible boundary 
conditions at the junctions are not yet fully understood. 
The dynamics away from equilibrium needs further inves- 
tigation as well. In this context we focus below on three 
particular aspects. These are the interplay between bulk 
dynamics and vertex interactions, the construction of an 
equivalent total scattering matrix relative to all external 
leads of the network and finally, the impact of an am- 
bient magnetic field on the conductance and the noise 
away from equilibrium. We consider the first problem in 
general and, keeping the discussion as simple as possible, 
illustrate the last two on the example of a ring junction 
in the presence of a magnetic flux (p, as shown in Fig. 
[2] More precisely, we investigate the conductance and 
the noise of this junction as a function of the transmis- 
sion t in the ring, the flux 0, the temperatures Pi and 
chemical potentials fii of the external edges, as shown in 
Fig. [3] below. Aharonov-Bohm type oscillations with 
occur in both the conductance and the noise. The period 
of these oscillations equals the elementary flux quantum 
(j)o — 2'Khc/e associated with a single charge e. The pure 
thermal {fii = 0) noise has a 0-dependent power law be- 



haviour at small temperatures, with discontinuities in the 
points (/) — n(j)Q , n G Z. As functions of /Xj, the current 
and the shot noise have an interesting plateaux structure 
in the regime t ^ 1. The fundamental and essentially 




FIG. 1: (Color online) A graph F with vertices Vi and edges 
Ej modeling the junctions and the leads of a quantum wire 
network. 
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FIG. 2: (Color online) Ring junction with local S-matrices 
{Si : 1 = 1,2,3} and magnetic flux (p (left) and its effective Y- 
junction counterpart with equivalent total scattering matrix 
(right). 



unique input for deriving these results is the requirement 
of self-adjointness of the Schrodinger and Dirac Hamil- 
tonians on the ring graph of Fig. [2] in external magnetic 
field. This requirement determines a family of admis- 
sible point-like interactions at the vertices on the ring. 
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We focus on the scale invariant ones because they incor- 
porate the universal features of the system while being 
simple enough to be analyzed explicitly. In particular, 
no phenomenological potentials are introduced. 




FIG. 3: (Color online) Y-junction connected at infinity to 
thermal reservoirs with inverse temperature /3i and chemical 
potential fj,i. 



The paper is organized as follows. In the next sec- 
tion we discuss the relation between bulk dynamics and 
point-like interactions at the vertices. The main tool here 
is the use of the self-adjoint extensions of the Hamilto- 
nian from the bulk to the whole graph. In section III we 
reconstruct the global scattering matrix S*^ from the local 
5i-niatrices in the presence of magnetic flux (p crossing 
the ring junction in Fig. [2] We describe the basic fea- 
tures of S*^ and establish the breaking of time reversal 
invariance induced by (j). In section IV we investigate the 
non-equilibrium dynamics of the Y-junction with exter- 
nal leads connected to thermal reservoirs with a differ- 
ent temperature and chemical potential, as shown in Fig. 
|3j The conductance and the noise are explicitly derived. 
The helicity transport in the Dirac case is also described. 
Section V is devoted to the conclusions and a discussion 
of some further developments. The appendices collect 
some technical details. 

Let us mention finally that the wires, displayed in the 
figures of this paper, are planar and in most of the cases 
have straight line edges. However, the discussion below is 
completely general and applies to segments of arbitrary 
smooth curved^ in M.^ as well. It is also extremely impor- 
tant to realize that in the mapping, presented in Fig. [2] 
the Y-junction with scattering matrix S"^ is a convenient 
way to model the ring junction with matrices Si and is 
not the result of taking a zero size limit of the ring while 
maintaining a fixed magnetic flux. In spite of the fact 
that we adopt scale invariant local scattering matrices 
Si, the total scattering matrix S"^ encodes explicitly the 
finite size of the ring and is not scale invariant. The ring 
junction in the left hand side of Fig. [2] and the point like 
Y-junction in the right hand side are completely equiva- 
lent and we use these names interchangeably in the rest 
of the paper with this understanding. In this re spect 
our approach is very different from the on^^SlMIlIl using 
boundary conformal field theory. 



II. INTERPLAY BETWEEN BULK DYNAMICS 
AND VERTEX INTERACTIONS 

The fundamental inputs for applying transport theory 
on graphs are essentially two: 

(i) dynamics in the bulk (edges); 

(ii) interaction at the vertices. 

In most of the papers on the subject it is implicitly as- 
sumed that (i) and (ii) are two independent ingredients. 
Some recent development^^USI] the spectral theory 
of operators on graphs ("quantum graphs") have shown 
however that this is not the case, if one requires unitary 
time- evolution of the system. It turns out that this basic 
physical condition relates the bulk and vertex interac- 
tions. The point is that the bulk dynamics is described 
by a Hermitian Hamiltonian, which becomes self-adjoint 
only by imposing special boundary conditions at the ver- 
tices. These conditions generate particular point-like in- 
teractions, which are described by specific (and not arbi- 
trary) scattering matrices Si, associated with each vertex 
Vi of the graph. Since the form of Si has relevant physical 
implications, it is instructive to describe the phenomenon 
in more details. We consider for this purpose the simplest 
example of fermion excitations, which propagate freely in 
the bulk. 




FIG. 4: (Color online) A vertex V with n edges Ei. 

Let us denote by a generic vertex of the network 
with n edges Ei, orientecP^l as in Fig. [ij We consider 
first non-relativistic fermions (Schrodinger network). In 
this case the time evolution in the bulk is described by 
the Schrodinger equatiorP^l 

{idt + :^dl^i;{t,x,i)^0, (1) 

where a; > is the distance from the vertex V. The only 
non-trivial interaction is localized at V and is encoded in 
boundary conditions at x = 0. To select these conditions, 
we require that the bulk Hamiltonian —9^ admits a self- 
adjoint extension at a; = 0. In such a way we cover all 
point-like interactions localized at V and leading to a 
unitary time evolution of the system. The most general 
boundary condition, implementing this natural physical 
condition for the Schrodinger network, iPSHiU 

n 

J2 - c/)^, ^(t,o, j) - *(i + u),j{d,m, 0, j)] = , 

(2) 



3 



where U is an arbitrary n x n unitary matrix and A is 
a real parameter with the dimension of mass. U = I 
and U = —I generahze to the star graph in Fig. |4] the 
famihar Neumann and Dirichlet boundary conditions on 
the half line. It has been established'^^ ''^^ that the point- 
like interaction, induced by generates 



S{k) = - 



[X{I - U) - k{I + U)] 
[X{I - U) + k{I + U)] 



k e 



(3) 



which is a very special momentum-dependent scattering 
matrix. An equivalent form of ([S]), which is more conve- 
nient for explicit computations, is given in the appendix 
[A| The representation (A3) shows in particular that S{k) 
is a meromorphic function with simple poles, all of which 
located on the imaginary axis and different from 0. 

The equation of motion ([T]) is invariant under the time 
reversal operation 



Tip{t, X, i)T 



-VT'iPi-t,x,i) , 



\Vt\ 



(4) 



T being an anti-unitary operator. We stress however that 
the vertex interaction (and therefore S{k)) preserves the 
time reversal symmetry if and only if U is symmetric!^, 
or equivalently 



S{kY = S{k) 



(5) 



where the apex t indicates transposition. We will apply 
the condition ([5| when discussing below the breaking of 
time reversal symmetry in the Y-junction with magnetic 
flux. 

Let us consider now relativistic fermions (Dirac net- 
work). In this case the dynamics is governed by the Dirac 
equation (we consider for simplicity the massless case) 



iltdt - lxdx)i^{t, x, i) = , 



x> 0, 



where 



and 



■)p{t,x,i) 



_ ( V'i(^,a;,i) 



It 



(6) 



(7) 



(8) 



The boundary conditions, which define all self-adjoint 
extensions of th e bulk Hamiltonian ijtjxdx at the vertex 
V, are noWSHU 



?/;i(i,0,i) =^C^yV'2(i,0,j), 



(9) 



where U is any unitary nxn matrix. Observing that both 
the equation of motion ^ and the boundary condition 
([9]) preserve scale invariance, it is not surprising that the 
scattering matrix, associated with (jol, is simpljB^ 



s{k) = e{k)u + e{-k)u-\ 



(10) 



9 being the Heaviside step function. The necessary and 
sufficient condition for time reversal invariance in the 
Dirac case is ([s]) as well. 

At this point we would like to draw the reader's atten- 
tion to a peculiar feature of the above systems, related 
to the fact that the argument i in the field ip{t, x, i) is a 
space coordinate and does not label an internal degree of 
freedom. For this reason the transformation 



S{k) 



S'{k)=US{k)U- 



(11) 



tp(t,x,i) I — > ^p'{t,x,i) = ^U^jil;{t,x,j) , (12) 

where U is a n x n unitary matrix, is allowed only for 
diagonal U. Otherwise, the field ^' in ( 12 ) is delocal- 



ized in different edges and therefore does not describe 
conventional physical excitations. 

Summarizing, we have classified and described explic- 
itly all point-like interactions at the vertices which lead 
to unitary time evolution of the Schrodinger and Dirac 
networks. Due to the non-trivial four-fermion bulk in- 
teractions, the case of a Luttinger network is more com- 
plicated. The possibilitjI^nHSIl to impose there non-linear 
vertex boundary conditions, formulated directly in terms 
of the electric current, is very attractive but needs further 
understanding and will be investigated elsewhere. 



III. S-MATRIX FOR THE Y-JUNCTION WITH 
MAGNETIC FLUX 

We consider below the ring junction (Fig. [2]) in the am- 
bient space M'^ with a time-independent electromagnetic 
field, generated by the potential A(x). The magnetic 
fiux, crossing the junction, is given by 



£a(x) 



(13) 



where the contour of integration C is the ring in Fig. 
[2] We indicate now the basic steps for reconstructing the 
total scattering matrix S"^ from the three local S'-matrices 
at the vertices on the ring. Several frameworks [IS1I5T| 
have been proposed for this purpose. We follow below 
the approach of Ref. Ell which is better suited to the 
case with ambient magnetic field and provides explicit 
expressions. The starting point is to implement in ([T]) 
and ^ the minimal coupling 



idx + eAxix, i) . 



(14) 



where Ax{x,i) is the projection of the vector potential 
A(x) along the edge Ei of the graph and e is the electric 
charge. For simplicity we set e — 1 throughout the paper. 
With this convention, the minimal coupling gives rise to 
(/)-dependent phases that particles pick up as th ey tra vel 
along the edges. These phases can be transferrecP^'^ by 
a kind of "gauge" transformation (affecting the boundary 
conditions (l2|l9|) to the matrix U and therefore to Si{k). 
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Since the explicit form of S'^ for a generic ring junction 
with general Si{k) is quite complicated, we simplify the 
considerations by focusing on the case of identical local 
scattering matrices Si{k) = S{k) (i = 1, 2, 3) and equidis- 
tant vertices, separated by a distance d along the ring. 
In this case the system is invariant under cyclic permu- 
tations, implying that S"^ is a circulant matrix, i.e. 

/ cri(fc,(/)) (T2{k,(p) a3{k,(f))\ 
§"*(fc)= a3(fc,0) ai(fc,0) (T2(fc,<^) . (15) 

yCT2(fc,0) 0-3(^,0) CTi(fc, (/>) J 



The explicit form of (Ji{k,(j)) in ter ms of the matrix ele- 
ments Sij{k) of the local S'-matrix i^SHMl 



,ik)^l±e^-i^~^^0~^^x(k,tt^^^) , (16) 



where 



A(fc,0) = -det[5(fc)] 



\kd 



^23(-fc)e'^ - S32i-k)e-'^ - dct[5(-fc)]5ii(fc). 
-''"^ - <S23(fc)e'<^ - S'32(fc)e-''* - dct[S'(fc)]S'ii(-fc)e"='' 



(17) 



For analyzing the universal features of the Y-junction, 
it is enough to concentrate on scale-invariant (critical) 
local 5-matrices. These matrices are classified in Ap- 
pendix |A] and, for the Schrodinger junction, form the 
two one-parameter families Sj. defined by equation ( A8 1 . 
Plugging (A8) in (17), one obtains 



except at the points 



^ ,\ t(cos(/)=F cosfcd) ± i(t — 1) sinfcd 

' ^' ~ ^^(cos (j) T cos kd) =F i(t - 1) sin kd ' ^ ' 

Equations ([Tsj [16] [Ts] ) determine the set 

{§,'i{k;d,t) : d>0,te [0,1]} (19) 

of total scattering matrices, parametrized by the size 
d of the internal edges and the transmission coefficient 
t among them. The boundary of the parameter space 
needs special attention. §^ simplify at d = 



St(fc;0,t) = Tl, 



(20) 



which describe three disconnected edges. For this reason 
we take d 7^ in what follows. The behaviour at i = 
0, 1 is responsible for a kind of resonant tunneling effect, 
discussed later on in this paper. One has in fact 



Si(fc;d,l) = Tl, S^(/c;d,0) =±I. 



(21) 



for generic values of fc, except at the values (see eq. (33)) 
corresponding to resonances described below. Moreover, 
one easily verifies that 



§^(A:;d,t) = -SV^+"^(fc;d,i). 



(22) 



Therefore, without loss of generality we concentrate in 
what follows on S^l, omitting for simplicity the index 

Let us discuss at this point in more details the depen- 
dence of S"^ on the magnetic flux and the momentum 
k. First of all we observe that 



— 3n7r , 



(24) 



Taking into account the condition ([s]), we thus conclude 
that only the fluxes ( 24 ) preserve time reversal invari- 



ance. Otherwise, the magnetic field breaks down the time 
reversal symmetry. In fact, the following relations holcP^l 



iS^)\k) = S'^(-fc) , (§^)*(A:) = §~^(fc) , 
so that 

a*{k,(j)) = (Tj(-fc, -0) , 

and 

cri(fc, (f>) = cri(fc, -<p) , (T3(/c, (j)) = a2{k, -(f>) 
We can now define the transmission amplitudes 

T+{k,(f>) = \a2{k,(l))\ , T^{k,(f>) = |cr3(fc»| 
and the reflection amplitude 

Q{k,(j)) = |cri(fc,0)| , 



associated with (15). They satisfy 

T±(fc,0) = T±(-fc, -(/)) , T_(/C,(/)) = r+(fc, 

and 

£i(fc, (j)) = g{~k, 0) = g{k, ~(f)) , 



(25) 
(26) 
(27) 
(28) 
(29) 
), 

(30) 
(31) 



which refiects the symmetry structure of our ring junc- 
tion. 

From the unitarity of S"^ we deduce the expected rela- 
tion 



(^2{k,(t>) ^ (T3{k, 



(23) 



g\k,cb)+Tl{k,^)+Tl{k,(l))^l. 
among probabilities. 



(32) 



5 



One can verify that both t± and g are periodic in (j) 
with period <^o = 27r. Note, that this is not the case 
of (72 and (T3, which are only 67r-periodic. The period 
00 has a deep physical meaning. In fact, recalling our 
convention c — H = e = 1, (pa equals precisely the el- 
ementary flux quantum 2'Khc/e associated with a single 
charge e. This flux quantum appears^*^ naturally in the 
context of Aharonov-Bohm type oscillations in non sim- 
ply connected mesoscopic systems. It is not surprising 
therefore that the 0o-periodicity shows up below in the 
currents, the conductance and noise of the ring junction. 

We emphasize finally that both g and t± are k- 
dependent, in spite of the fact that the local scattering 
matrices (A8) are constant. This dependence is a di- 
rect consequence of the finite size of the ring. In fact, 
the momentum k enters g and t± only through the di- 
mensionless combination kd. It follows from (17 1 that 



g and t± are 27r/c?-periodic in fc, the shape of the oscil- 
lations being strongly influenced by the transmission t 
and the flux 0. The behaviour of the probability for 
d = 1, shown in Fig. [5] confirms this statement. The 
dashed (black) and the continuous (red) lines describe 
the oscillations for t = 0.5 and t = 0.99 respectively. As 
already observed, for i ~ 1 the external edges are almost 
isolated. Accordingly, in this regime, one expects very 
small transmission amplitudes. This is indeed the case 
with the exception of the momenta k — ^'^^J'^" , charac- 
terized by the appearance of sharp peaks with maximum 
close to 4/9 and corrections of order 1 — t. For (j> ^ titt, 
there are six of them in each k interval of length 27r/c? 
and only three ii 4> — mr. Fig. [5] illustrates the phe- 
nomenon for = (left) and = 7r/4 (right). Because 
of (32 1, the behaviour of the reflection amplitude g is 



complementary. 
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FIG. 5: (Color online) ri(fe,0,d = 1) for t = 0.5 (dashed 
black line) and t — 0.99 (continuous red line) for (j) ~ (left) 
and (j) = n/4: (right). 

We stress that at t = 1 the amplitudes t± become 
actually discontinuous in k. One has indeed 



k = ^^^, nez. 




The impact of the peculiar A:-dependence of t± for t ^ 1 
on the current and the shot noise is discussed in the next 
section. 

It is worth mentioning that similar phenomena show up 
in the other limit t ^ of almost disconnected external 



edges. Also in that case the transmission amplitudes are 
negligible, apart around the points k = nn, where the 
nontrivial behaviour, displayed in Fig. [6] is observed. 
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FIG. 6: (Color online) rl{k,(t),d = 1) for t = 0.01 around 
k — (left) and k — tt (right). The dashed (black) and 
continuous (red) curves correspond to (ji = and (p — n/A 
respectively. 

Another type of discontinuities of t±, which involves 
the magnetic flux and also deserves attention, is de- 
scribed by 



lim T± (fc, (j)) 



I , = n(f>Q . 
, 4>^ n(f>o . 



I, n & Z , 



(34) 



We will show below that ( 34 1 has important consequences 



for the behaviour of the pure thermal noise at small tem- 
peratures. 

Let us conclude this section with a remark about the 
Dirac junction. According to the observation at the 
end of Appendix A the matrices §^ apply both to the 
Schrodinger and Dirac junctions, the possible differences 
in the behaviour of the conductance and the noise in these 
cases being related essentially to the different dispersion 
relations. 



IV. CURRENTS, CONDUCTANCE AND NOISE 
AWAY FROM EQUILIBRIUM 

In this section we apply the results of Ref. [33] about 
non-equilibrium transport on star graphs to the Y- 
junction in Fig. [3j The interaction in the vertex is de- 
scribed by the flux dependent scattering matrix and 
the leads are connected at infinity to thermal reservoirs 
(heat baths) with (inverse) temperature /3i and chemical 
potential 



fJ-i 



(35) 



where kp defines the Fermi energy and Vi is the ex- 
ternal voltage applied to the edge Ei. In what follows 
we keep kp fixed, varying eventually the gate voltages 
Vi. The system is away from equilibrium if admits 
at least one non-trivial transmission coefficient among 
edges with different (Hi and/or /i^. It turns outP^l that 
the corresponding non-equilibrium dynamics can be im- 
plemented by a steady state fJ^g.^, characterized by non- 
vanishing time-independent charge and heat currents cir- 
culating along the leads. Nevertheless, the contact with 
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the heat baths ensures energy conservation. In this sense 
describes an ideal system without dissipation. The 
construction of i^^,^ involve^^the scattering matrix S"^ 
and fully takes into account both the minimal coupling 
with the external magnetic field and the vertex inter- 
actions. We denote in what follows the expectation 
values in the state il^,^ by {■■■)p^fi and consider the 
Schrodinger and the Dirac cases separately. We stress 
that the current correlation functions {jx{t,x,i)) p^^j, and 
{jx{ti,Xi,ii)jxit2,X2,i2))i3,ii used below are exact. No 
approximations, like linear response theory, are adopted. 



A. Schrodinger junction 

1 . Current 

The charge transport in this case is associated with the 
conserved current 

jx{t,x,i) = ^ [ipHdx'4') - {dxi^'')i^] {t,x,i) . (36) 



Using the explicit constructiorP ^of th e steady state ri/j,;^, 
one finds the Landauer-Biittikei'^^'^ formula 



dk k 
27r m 



where 



g-/3,[c^(fe)-Aii] ^2 



(37) 



(38) 



is the familiar Fermi distribution. Unitarity of S"^ implies 
the Kirchhoff 's rule 



^{jx{t,x,i))p^i, = 0. 



(39) 



Introducing the variable ^ = k'^ /2m and using that the 
Fermi distribution (38) approaches the Heaviside step 



bmction 9{pLi — ^) in the zero temperature limit /?; — cx), 
one obtains from (|37l 



Ji = {3x{t,x,i))^ 

3 



^(M.)g-E^(^.)r^K(v^)r- (40) 



We will investigate the behaviour of (40), comparing J, 
to the shot noise in subsection 3 below. 



and (37) determine the conductance tensor 
I n dk 



mVi 



3 Jo 



27r 



\St^{k)\' d,{k). (42) 



The periodicity of |cri(fc)| in 0, established in the previous 
section, indicates that Qij oscillate with period (pQ. The 
unitarity of S'^ implies the following ^-independent bound 



dk 
2tt 



kdj{k) 



1 



log 



2nP,\V,\ 

on the amplitude of the oscillations. 



(43) 
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FIG. 7: (Color online) Gii(<^,t) for /3i = 1 (left) and /3i = 0.2 
(right). 
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FIG. 8: (Color online) Sections of the plots in Fig. [t] 
4> (first line) and fixed t (second line) at /?i = 1(1 
/3i = 0.2 (right). 



for fixed 
left) and 



2. Conductance 



The standard definition 



In order to get a more precise idea on the dependence 
of Gi, on (j), the transmission t and the temperature /3, 



we concentrate on (42). The fc-integration can not be 



performed in a closed analytic form, but being well de- 
fined, the integral can be computed numerically. The 
(41) following plots illustrate the result of this computation, 
in which the values d = 1, m = 1/2, kp = 2 have been 
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FIG. 9: (Color online) Gi2(0, at j3i — 0.2 and its sections 
at t = 0.3 (dashed line) and t — 0.8 (continuous line). 

fixed. Typical contour plots of Gn in the (j) — t plane are 
displayed in Fig. [7]at two different temperatures (/3i = I 
and /?! = 0.2) for the voltage ^1 = 1. As usual, higher 
regions are shown in lighter shades. Fig. [8] collects some 
sections of the contour plots. The first line displays the 
sections at </> = (dashed line) and = 27r/3 (continuous 
line). As expected, the conductance vanishes at t = 0,1, 
when the external edges of the Y-junction are isolated 



from each other. We see also that the position of the 
maximum of Gn depends on the flux 0. In the second 
line of Fig. Islwe report the sections at t = 0.4 (dashed) 
and t = 0.7 (continuous), which show the expected os- 
cillation in with period (j)o. The continuous (red) lines 
illustrate the impact of the higher harmonics n(j)Q, which 
become relevant for t > 0.5. Note that the oscillations of 
Gn are symmetric with respect to (/< = tt, which is actu- 
ally a general feature of the reflection amplitude g. The 
off-diagonal elements do not share this property, which 
is evident from Fig. [9| showing the behaviour of G12. 

3. Noise 

The zero frequency noise power is defined as usuaf^^ 

by 

/oo 
die-* {j,{t, xi, z)j,(0, a;2, j))^°7 , (44) 
-00 

where {jxiti,xi,i)jx{t2,X2, is the connected 

current-current correlation function in the state J^^,^. It 
turns ouiP^that Py is a;i^2-independent and is given by 



m 



1 c\k r 

ViJo ^''{^^^d,(k)c,{k) - \St^{k)\^d,{k)c,{k) - \S^^{k)\^d,{k)c.{k) 
\ ^ Sf,(fc)§5(fc)§|„(fc)st(fc)h(fc)rf™W+c„(fc)d,(fc)]}, (45) 



/,m— 1 



where Ci(fc) = I — diik). P is a symmetric matrijP^J. If 
we assume now fjn = ^ and /3i = /3, so that the setup of 
Fig. [3] respects the cyclic symmetry among the vertices 
on the ring, then P is also a circulant matrix. Combining 
this fact with the Kirchhoff rule 



we obtain that 



-0, 


i — 


1,2,3, 


(46) 


I 


-1/2 


-1/2 \ 




1/2 


I 


-1/2 , 


(47) 


1/2 


-1/2 







with 



Pu = 



m 



dk 
2t: 



kd{k)c{k) [Tl{k,^)+Tl{k,(^)\ 



(48) 

Like the conductance, Pn oscillates in with period ^g. 
The bound on the amplitude, following from unitarity is 



now 



< Pii < 



Ah. 

^^kdik)cik)^-^ 



m 

Let us investigate the pure thermal noise, obtained 
from (48 1 by setting /i = (i.e. a voltage V — kp). Using 
the variable T ~ 1//?, the numerical study confirms that 
at large temperatures 

Pn-T, for T^oo, (50) 

independently of the flux cj), as suggested by ( 49 1 . In this 
regime one recovers therefore the well-known Johnson- 
Nyquist formula. The situation changes drastically at 
small temperatures, where we flnd a 0-dependent power 
law behaviour 



10 



for P^O, (51) 
In absence of magnetic 



where g is plotted in Fig. 
flux one observes a linear behaviour at small tempera- 
tures as well. However, as soon as the magnetic flux is 
switched on, the exponent g{(f>) rapidly increases from 



8 



3 4 



FIG. 10: (Color online) The power g{(j>) in the Schrodinger 
case. 



/? — CX3 limit, keeping fii > arbitrary. Adopting the 
variable ^ = k'^ /2m one gets 



where 



J-,,,™(fc) = Sf,(fc)sJ(fc)S^,„(fc)sL(A;) , (57) 



and £ is the sign function. Note that since Pij is sym- 
metric and satisfies ( |46[ ), we only need to compute the 
diagonal elements {Pa : i = 1,2,3} in order to recon- 
struct the complete matrix. Indeed, for i, j, k mutually 
distinct, 



the value 1 to 2, leading to a quadratic power low. We 
see therefore that the presence of a magnetic flux in the 
Y-junction implies a relevant small temperature modi- 
fication of the standard Johnson-Nyquist thermal noise 
behaviour. This new feature provides an interesting sig- 
nature of a physical effect that hopefully can be observed 
experimentally. 

The above numerical results are supported by an exact 
analytic argument, reported in a more formal accompa- 
nying paper 



which gives 



7^ rKpo . 



neZ, 



(52) 



(58) 



The diagonal elements read 

3 



Pa = ^ £(A*; - ^m) / ^J^iilrn (^\/2m^j . (59) 



/<m=l 



Assuming for definiteness that /ii < /i2 < /is, one obtains 
from (l59| 



Pu 



d^ 

, 2^ 



2tt 



tI{1-tI), (60) 



Strictly speaking, the power g is discontinuous in the 
points (j) — n(f>Q. One can sho'vJ'^ that the origin of this 
discontinuity is the peculiar behaviour of the transmis- 



sion probabilities in the integrand ( 48 1 , namely the 
lack of commutativity of the two limits 



lim lim Tj.{k) ^ lim lim 



r|(fc). 



In fact, it follows from (34) that 



lim 



0— ^n0o k^O 



lim r^(fc) = 0, 



whereas 



lim lim t. 



l(fc) = 



(53) 



(54) 



(55) 



The apparently smooth behaviour of g{(t>) around the 
points (p = 0,(pQ in Fig. [TO] is clearly an artifact of the 
numerical computation, which is performed unavoidably 
with a finite precision. 

The situation is more complicated in the presence of 
a non-zero chemical potential /i > 0. The parameter 
controlling the behaviour of the noise Pu in this case is 
n/T. For the complete study we refer to Ref. [59], where 
various ranges of this parameter are considered. 

Let us investigate now the shot noise defined by 
eq.(45). For this purpose we set /3i = (3 and take the 



P22 — 



-4), (61) 



P 



33 



2tt 



<(l-<) + 



^9'il-9'), (62) 



where g and t± are computed at fc = \/2m^. The be- 
haviour of g and t± implies that Pa oscillate in cf) with 
period The amplitude is subject to the obvious uni- 
tarity bound 



0< Pti< H3- Ml 



(63) 



The plots of the shot noise Pu in Fig. [TTJ obtained 
for /ii — I, fi2 = 2, fi3 — 3 and m — 1/2 confirm the 
oscillations and the impact of the higher harmonics at 
t ~ 0.8 (continuous red line). 

We study finally the behaviour of the shot noise as a 
function of the chemical potentials /li, or equivalently, the 
voltages Vi in (35). It is instructive to do this, comparing 
Pu with the zero-temperature steady current Ji given 
by (40), and the transmission amplitude t±. For this 



purpose we fix m = 1/2, ni = fi2 = d = 1 and vary /is. 
In this regime 



27r' 



J3 = - Ji - J2 , (64) 
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the next in Fig. 12 for = 7r/4, in correspondence with 
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FIG. 11: (Color online) Pii{(j),t) and its sections at t — 0.3 
(dashed line) and t — 0.8 (continuous line). 



the fourth peak in r_. However, if </) < 7r/4 roughly, the 
peak is moved to the left, so that Ji takes on the higher 
value. If (/) > 7r/4 roughly, the peak moves to the right 
and Ji takes on the lower value. 
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FIG. 13: (Color online) Plot of Ji at t = 0.99 and ^3 = 20 as 
a function of 6. 



B. Dirac junction 



FIG. 12: (Color online) Plots of Pn (black dashed), Ji (blue 
dotted) and rLjl-n (red continuous) at t = 0.99 for = 
(left) and = 7r/4 (right). 



'^^r!(l-r2) = P22M), (65) 



2tt 

„2/ 



(66) 



The Dirac junction (|6| |9j) can be treated in the same 
way as the Schrodingcr one. For this reason we omit the 
details and report only the basic equations, discussing 
the differences. First of all, in the Dirac case both parti- 
cles (with chemical potential fii) and antiparticles (with 
chemical potential Jli) are present. The relative steady 
state i^i3,fj.,ji has been constructed in Ref. ESI The charge 
transport is associated with the current 



{t, X, i) = [. ^l^i : - : tljlil)2 :] (i, x,i), (67) 



An interesting resonant tunneling effect is observed for 
t ^ \. This corresponds to the situation where the exter- 
nal edges are weakly coupled to the ring. The peaks in 
the transmission amplitudes r-t, discovered in section III, 
can be interpreted as resonances corresponding to eigen- 
states of the ring. A similar situation was discussed irP^ 
in the case of the ring with two external edges. As the 
voltage is increased, these resonances generate plateaux 
in the shot noise Pa and the current Ji. This fact is il- 
lustrated in Fig. 12 where we plotted t^{^/JI^) /2tt (con- 
tinuous red curve) , Pii{^^) (dashed black curve) and 
Ji(/i3) (dotted blue curve). Switching on the magnetic 
field changes the location of the peaks and hence the lo- 
cation of the jumps from one plateau to the next. 

For a fixed potential /.ta, the amount of current flowing 
in the external edges can vary with the magnetic flux. 
This is shown in Fig. [13] We fixed /i3 = 20 which is the 
approximate value where Ji jumps from one plateau to 



where : • • • : indicates the normal product. Introducing 
the Fermi distributions for particles and antiparticles 



-/3.(|fe|-M.) 



~ft(lfcl+/I.) 



1 + e-ft(|fe|+M.) ' 
(68) 



one finda^ in this case the following steady current 







dfc 

2^^2-1^^^ 



[%~|Sf/A:)p] [/,(fc)-/,(fc)], (69) 



conductance 



°° dk r 

2Tr 



sf,(fc)p im-m], (70) 



and zero frequency noise power 
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g{%F,,(fc)-|Sf^.(fc)pF„-(fc)-|S^,(fc)P^^,,(fc) + i ^tlimU''PU>')^trnimFUk)+Frnlik)]} , (71) 



l.m—1 



r 



where 

F^,{k) = Mk)[l - f,{k)] + Mk)[l - f,{k)] . (72) 



The heUcity transport, which has attracted some atten- 
tion recentlji^, is described in our case by 



+ :4i:2-]{t,x,i). (75) 




Note the change of sign of the second term in the right 
hand side of ( 75 ) with respect to the charge current ( 67 1 . 



This sign difference propagates in the hehcity steady cur- 
rent, which reads 



dfc 



|§f/fc)p] [/,(fc)-/,(A:)] 



(76) 



Unitarity imphes 



3 1^1 

FIG. 14: (Color online) The power h{(j)) in the Dirac case. ^~^(fe^(t, x, i)) p^fx.jj — — — 



vr — ft 



In 



1 



(77) 



Apart from the antiparticles and the hnear dispersion 
relation in the Fermi distribution, the only difference be- 
tween ( |70|71| ) and ( |42|45[ ) is a factor k/m m the inte- 
grand. For this reason, the behaviour of the Dirac junc- 
tion is similar to that of the Schrodinger one. One excep- 
tion of this rule is the small temperature behaviour of the 
shot noise Pn, following from (71 1 by setting /i^ = /ii = 
and Pi = j3. A computation, analogous to that in section 
IV. 3, gives now 



which shows that the relative Kirchhoff rule is not auto- 
matically satisfiecP^. Differently from the electric charge, 
which is always conserved for a unitary S"^, the helicity is 
a conserved quantum number in the state if and 

only if the temperatures and the chemical potentials of 
the heat baths satisfy 



1 



1 



0. 



(78) 




with 



for 
for 



4> = "00 : 



T 
T 



oo , 
0, 



n e z. 



(73) 



(74) 



The result of the numerical computation, plotted in Fig. 
[T4l confirms this statement. 

We see that, as in the Schrodinger case, magnetic fluxes 
(j) ^ n(j)Q modify for T ^ the usual linear temperature 
dependence of the pure thermal noise. In spite of the 
fact that the scattering matrix S"^ for both Schrodinger 
and Dirac junctions is the same, the difference of the 
dispersion relations is clearly visible from the critical ex- 
ponents g{(j)) and h{(j))^ which for (p ^ rKpQ equal 2 and 3 
respectively. In principle, this fact could be used exper- 
imentally to determine whether a given junction is best 
described by the Schrodinger or the Dirac equation. 

Besides a conserved electric current, in the massless 
Dirac case one has also a conserved helicity curreni!^. 



1 + e-i^-t^ 

This is for instance the case if fii = ~Jli with generic ft 
V. OUTLOOK AND CONCLUSIONS 



Combining a scattering approach of quantum field the- 
ory with results arising from the spectral theory of op- 
erators on graphs, we discussed transport properties of 
fermions on a ring connected to thermal reservoirs. In the 
bulk the fermions interact only with an external ambient 
electromagnetic field. At the vertices on the ring we con- 
sider the most general scale invariant local interactions 
which generate unitary time evolution. The formalism 
applies equally well to the Schrodinger or Dirac case and 
conceptually uses only the self-adjointness of the relative 
Hamiltonians on the ring graph in Fig. [2] with magnetic 
flux (/). 

As is well-known, the crucial ingredient in the scatter- 
ing approach to transport theory is the scattering matrix 
of the "sample" under investigation. In our case, this 
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matrix takes into account all the above mentioned inter- 
actions. The systematic method, originally developed in 
[5T] , can be generalizecP^ to the presence of an external 
magnetic field and provides us with a method to compute 
the relative scattering matrix S"^ explicitly. Hence, all the 
required analytic structure of reflection and transmission 
amplitudes is known and its dependence on and the 
size of the ring is easily incorporated. 

With this last ingredient now available, this paper rep- 
resents a further step in exploring the interesting in- 
terplay between spectral theory on graphs and quan- 
tum field theory on graphs with applications to con- 
densed matter problems, as initiated in the series of 
papers ^3(1^^33j. The non equilibrium dynamics, gener- 
ated by the contact of heat baths with different tem- 
peratures and chemical potentials, is captured by steady 
states ^i3,fi, which incorporation ^jjg scattering matrix 
S''^. The expectation values of the current in fi^g,^ re- 
produce thfiLandauer-Biittiker steady current and noise 
formulai^^MIII] ^gpj-^s of Using these formulas, we 
found a resonant tunneling effect in the case where the 
ring is weakly coupled to the external leads. The case 
with two external leads was originally discussed in [57] 
and we confirmed that the effect takes place also for three 
leads. In fact, using our general method, it is not hard to 
see that this effect will persist with more external leads 
[55] , the number of resonances available growing with the 
number of leads. It is worth mentioning that the reso- 
nant tunneling is a finite size effect. This is a reason why 
it has not been detected in the boundary conformal field 
theory approaclP^'^^'^ to Y-junctions. 

The most striking effect, discovered in this paper, is 
the sharp influence of the magnetic flux on the tempera- 
ture power law of the thermal noise and the typical sig- 
nature of the Schrodinger and Dirac cases. The usual 
Johnson-Nyquist linear behaviour is radically modifled 
at small temperatures to a quadratic (Schrodinger case) 
or cubic (Dirac case) behaviour as soon as the magnetic 
flux is switched on. That could not have been seen with 
previous methods available, since the effect is based on 
a detailed examination of the limiting behaviour of the 
fundamental eigenvalue of the scattering matrix S"^. 

As already explained, this paper represents a bench- 
mark where the tools of quantum field theory on graphs 
have been successfully tested and shown to produce re- 
liable and new results with relatively little effort. This 
opens the way to a range of future investigations. One 
is the study of quantum transport in more complicated 
three dimensional multiply connected networks. In prin- 
ciple, this can be done as easily as the present case, the 
essential features being captured by the scattering ma- 
trix of these structures. Another important extension 
is the use of bosonization techniques to describe quan- 
tum transport of self-interacting fermions in the bulk of 
the graph. This has already been implemented success- 
fully for the Tomonaga-Luttinger model on a star-graph 
[20121] . The next natural step is to investigate such a sys- 
tem on a ring. In this respect, it it important to realize 



that the present framework does not rely on conformal 
symmetry, thus allowing to investigate finite size effects, 
which is usually problematic with conformal techniques. 
A quantitative comparison with the exotic phase diagram 
of the Tomonaga-Luttinger liquid on graphs, studied in 
Refs. [10|22|24|27)28|30] . is an interesting issue for future 
investigation. 



Appendix A: Critical local ^-matrices 

Using the general expression ([s]), we derive here all 
scale-invariant (critical) 3x3 local S'-matrices for which 
the Y-junction in Fig. [2] is invariant under cyclic permu- 
tations. Let us denote by U the unitary matrix diagonal- 
izing U and parametrize 



as follows 



Ud = diag (e- 



2iai „2iQ2 



(Al) 



(A2) 



By means of ([3|, one easily verifies that U diagonalizes 
also §(fc) for any k and that 



Sd{k)=U^S{k)U = 

' k -\-irji k + i?72 



diag 



A; — iryi ' /c — i?72 ' 



where 



= Atan(ai) , 



2 - 2 



(A3) 



(A4) 



The scale invariant (critical) S'-matrices are k- 
independent and can be obtained by taking the limits 
rji or rji oo in (A3). Therefore, the most general 



scale-invariant scattering matrix, compatible with a uni- 
tary time evolution of the Schrodinger junction, is given 
by 

S^=USdU\ 5d = diag(±l,±l,.-.±l), (A5) 

where U is a generic n x n unitary matrix. From the 
group-theoretical point of view, any S'^ is a point in the 
orbit of some Sd under the adjoint action of the unitary 
group U{n). Obviously, one can enumerate the edges 
in such a way that the first p eigenvalues of S'' are — 1 
and the remaining n — p are -f-l. The cases p = and 
p = n correspond to S"^ = I„ and S"^ = — 1„ describing 
the Neumann and Dirichlet conditions respectively. 

The representation implies that besides being unitary, 
S'^ is also Hermitian, 



(A6) 



Assuming time-reversal symmetry in the individual local 
junction, according to ^ one has in addition that S"^ is 
symmetric. 



(5" 



5^ 



(A7) 
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Combining ( A6|A7 | with unitarity we conclude that the 
critical Schrodinger S'-matrices, preserving time-reversal 
invariance, are the symmetric matrices belonging to the 
orthogonal group 0{n). 

Let us concentrate finally on the case n — 3 (relevant 
for the Y-j unction in Fig. 2| and let us assume that the 
internal edges of each local junction are equivalent as far 
as transmission and reflection are concerned. Labeling 
these edges by the indices 2 and 3, the matrix elements 
{Sfj : i,j = 1,2,3} must be invariant under the ex- 
change 2 i — > 3. All these requirements fully determine 
two one-parameter families in 0(3), which can be repre- 
sented by 

/ l-2t ^2t(l - t) ^2t{l - t) 
S'|(t)=± y/2t{l - 1) t-l t 

\ y^2t{l - t) t t-l 

(A8) 

where t g [0, 1] is the transmission coefficient control- 
ling the local tunneling between the edges 2 and 3 of the 
junction. Since det{Sj.) = ±1, the matrices (A8) belong 
to the two disconnected components of 0(3) and to the 
orbits p = 2 and p = 1 respectively. 

The matrix Stit = 1/2) has been introduced in Ref. 
|56j . In Ref. [ST] the matrices Sj.{t) have been consid- 
ered for generic t g [0, 1]. To our knowledge, the above 
interpretation of Sj. (t) as critical points in the set of 
all scattering matrices, ensuring self-adjointness of the 
Schrodinger Hamiltonian on a star graph, is new. 

Finally, let us observe that Sj. = {Sj.) . Therefore, 
repla cing C/ in ( 10 ) by Sj. one concludes that the matrices 



( A8 ) describe critical points in the Dirac case as well. 



Appendix B: Comments on the construction of S** 

The main idea of the general method of Ref. "STto com- 
pute the total scattering matrix of an arbitrary quantum 



graph from the knowledge of its local scattering matrices 
and of its metric structure can be summarized schemat- 
ically as follows (for details in the case with an ambi- 
ent magnetic field see |55^). One collects the elements 
of the local scattering matrices in four matrices 5**^'""'"^ 
^(in-out)^ ^.(out-in) ^^^^ ^.(out-out) depending on whether 

the element connects two internal edges, an internal edge 
and an external edge or vice-versa, or two external edges. 
One also collects the information about the metric struc- 
ture of the graph through its consequences in terms of 
phases that the propagating modes pick as they travel 
along internal edges. In the presence of a magnetic field, 
one must also take into account the phase acquired by 
the propagating modes along the internal edges. All this 
information is gathered in a matrix E. The form of the 
total scattering matrix is then 

g ^(out — out) _|_ ^(out— in) ^(in— in)^ — 1 ^(in— out) 

(Bl) 

which has a very natural interpretation. An incoming 
mode on an external edge j is partly refiected straight 
after scattering off the vertex to which the external 
edges is attached (5*1°"*"°"'^) and partly after enter- 
ing inside the graph, undergoing all possible internal 
reflections and being emitted again on the same exter- 
nal edge ([S'(°"*-'")(i;-5('"-'"))-iS'('"-°"*)]^.p. The 
same mode is also partly transmitted to the ex- 
ternal edge £ after all possible internal reflections 

(|^5'(out-in)(i5 _ 5'(in-in)^-l^(in-out)j 

As expected, the symmetries of the graph allow for 
simpliflcations in the actual calculation of the total scat- 
tering matrix. This is shown explicitly in the case when 
the vertices form a regular polygon in |59| . In this paper, 
we deal with a particular case of this general situation 
where the three vertices form an equilateral triangle. As- 
suming identical local scattering matrices Si{k) = S{k) 

-J^that S"^ 



for alH = 1, 2, 3, one can s how" 
matrix given by ( 15|16|17). 
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